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ABSTRACT 

We report the discovery of two new low-mass, thermally bloated, hot white dwarfs among the Kepler 
sample of eclipsing binaries. These are KIC 9164561 and KIC 10727668 with orbital periods of 1.2670 
and 2.3058 days, respectively. The current primary in both systems is an A star of about 2 Mq. This 
brings the number of similar binaries among the Kepler sample to six, and the two new systems have 
the shortest orbital periods among them. The white dwarf in KIC 9164561 has the largest thermal 
bloating, compared to its cold degenerate radius, of about a factor of 14. We utilize RV measurements 
of the A star in KIC 9164561 to determine the white dwarf mass rather accurately: 0.197± 0.005 Mq. 

The mass of the white dwarf in KIC 10727668 is based on the Doppler boosting signal in the Kepler 
photometry, and is less accurately determined to be 0.266 ± 0.035 M©. Based on the inferred radii 
and effective temperatures of these two white dwarfs we are able to make an independent theoretical 
estimate of their masses to within ~0.01 Mq based on evolutionary models of their cooling history 
after they lose their hydrogen-rich envelopes. We also present evidence that there is a third body in 
the KIC 9164561 system with an orbital period of 8 — 14 years. 

Subjeet headings: techniques: photometric—stars:activity—binaries (including multiple)—binaries: 
general—stars: white dwarfs 


1. INTRODUCTION 

In relatively close primordial binaries which contain at 
least one star that can evolve substantially by the current 
epoch, mass transfer from the primary to the secondary 
is assured. The course of the evolution of each system de¬ 
pends on its initial orbital period, the two stellar masses, 
whether the binary is part of a triple system, spin and 
tidal interactions, and so forth. Some of the possible 
evolutionary channels have been explored and even sys¬ 
tematically mapped out by Rappaport et al. (1994), Iben 
& Tutukov (1995), Nelson & Eggleton (2001), Han et 
al. (2002; 2003), Podsiadlowski et al. (2003a), and Nele- 
mans (2010). Depending largely upon the mass ratio 
and the initial orbital period, the mass transfer may be 
stable and mostly conservative, stable and mostly non¬ 
conservative, or dynamically unstable. 

In the case of a close binary that comprises a primor¬ 
dial primary star with mass Mi^ < 2.5 Mq and a sec¬ 
ondary star that has M 2 ,o ^ 0.7 Mq, the mass transfer is 
most likely to be stable. However, depending on the ratio 
of thermal timescales of the two stars, the mass trans- 
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fer may range from highly non-conservative for small 
M 2 fl/Mifi ratios to mostly conservative for values of 
M^fl/Mifi close to unity. Whether the orbit ultimately 
evolves to be longer or shorter than the original orbital 
period depends on how much specific angular momentum 
is carried away with the matter ejected from the binary. 

In principle the best way to model such binary evolu¬ 
tion is via a stellar evolution code that can simultane¬ 
ously handle the nuclear and mass loss/gain evolution of 
both stars as well as the orbital evolution. The recently 
developed MESA code (Paxton et al. 2011) is capable of 
doing such binary evolution studies, but, to our knowl¬ 
edge, no systematic studies have yet been carried out of 
ordinary close primordial binaries. 

We know empirically, and from some theoretical stud¬ 
ies, that such evolution can lead to systems such as Algol 
binaries (see, e.g.. Nelson & Eggleton 2001) where the 
mass transfer is stable and the envelope of the primary 
has not yet been exhausted; sdB binaries (see, e.g., Han 
et al. 2002; 2003) which may have evolved from unstable 
mass transfer in a common envelope; systems with hot, 
bloated, low-mass white dwarfs (see, e.g., van Kerkwijk 
et al. 2010; Carter et al. 2011), where the envelope of the 
primary has been completely removed; and even double¬ 
degenerate binaries (Marsh et al. 2004; Iben, Tutukov, 
& Yungelson 1995; Nelemans 2010). 

In low-mass binaries where the mass transfer proceeds 
stably, regardless of whether it is or is not conservative, 
there is a fairly tight relation between the white dwarf 
mass, Mwd, and the final orbital period, Porb (see, e.g., 
Joss & Rappaport 1983; Savonije 1983; Paczynski 1983; 
Savonije 1987; Rappaport et al. 1995; Fig. 5 of Lin et 
al. 2011; Tauris & van den Heuvel 2014). For orbital 
periods 2 days < Porb ^ years, and white dwarf masses, 
A/wd 7 up to 1.4 Mq, in systems where the initial masses 
of the primaries, Mi^, were no larger than ^2.5 Mq, 
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the following expression gives the final orbital period in 
terms of the white dwarf mass: 


-Rorb 


1.3 X lO^Afe^^a 

(l + 4<d)''' 


days 


( 1 ) 


For orbital periods up to ^20 days, which apply to sys¬ 
tems with primaries of initial masses up to ~2.5M0, a 
more recently derived expression (Lin et al. 2011) is: 


.^orb 


4.6 X 10® 

(1 + 25M3| + 29M6 J3/2 


days . 


( 2 ) 


These expressions have been checked against a number 
of binary radio pulsars (Fig. 3 of Rappaport et al. 1995; 
Fig. 2 of Tauris & van den Heuvel 2014), and against a 
number of Kepler binaries containing low-mass hot white 
dwarfs (van Kerkwijk et al. 2010; Fig. 6 of Carter et 
al. 2011). 

Prior to this work, only four low-mass white dwarf sys¬ 
tems in the Kepler field that can be used to study this re¬ 
lation between orbital period and white dwarf mass were 
identified. They are KOI 74, KOI 81 (Rowe et al. 2010; 
van Kerkwijk et al. 2010), KHWD3 (KOI 1375; Carter 
et al. 2011), and KOI 1224 (Breton et al. 2012). In this 
work we report the discovery of new low-mass, highly 
thermally bloated white dwarfs in the Kepler binaries 
KIC 9164561 and KIC 10727668. The orbital periods of 
these binaries are 1.267 and 2.305 days, the shortest pe¬ 
riods among the six Kepler binaries containing low-mass 
hot white dwarfs. The white dwarf in KIC 9164561 is by 
far the most thermally bloated of any of the six white 
dwarfs. 

In Section [2] we discuss the search that led to these 
two systems among the Kepler binaries, and indicate how 
the data were prepared for analysis. Folded light curves 
are presented in Section along with a simple geometric 
analysis, including a harmonic decomposition of the out- 
of-eclipse portions of the light curve. The masses of the 
primary A stars are inferred from information on their 
temperatures and mean densities in Sect. In Sect. 
we present the radial velocity and spectral studies of KIC 
9164561. The mass ratio of KIC 10727668 is examined 
in the context of the Doppler boosting effect in Sect. 
Inferences about the mass ratios in the two systems are 
drawn from evaluating the ellipsoidal light variations in 
Sect. In Sect. |^the system parameters for these two 
binaries are evaluated with emphasis on the white dwarf 
properties. Evidence for a third star in KIC 9164561 
is presented in Sect. Finally, in Sect. [T^ we discuss 
KIC 9164561 and KIC 10727668 in the context of differ¬ 
ent binary evolution scenarios and white dwarf cooling 
models. 


2. DATA PREPARATION AND SEARCH 

In the process of searching for third-body companions 
to the binary stars in the Kepler sample (see, e.g., Rappa¬ 
port et al. 2013; see also Conroy et al. 2014), we produced 
both an O — C curve for the eclipses (“observed minus 
computed” eclipse times) as well as a folded light curve 
for each of the ^2600 cataloged Kepler binaries (Mati- 
jevic et al. 2012; Kirk et al. 2015)[^ We used the Ke¬ 
pler long-cadence time-series photometric data (29.4 min 
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samples) obtained from quarters 0 through 17. For all the 
known binary systems for which at least one quarter of 
photometry is available, the version 5.0 FITS Hies, which 
are available for all quarters, were downloaded from the 
STScI MAST websit43 We used the data that had 
been processed with the PDC-MAP algorithm (Stumpe 
et al. 2012; Smith et al. 2012), which is designed to re¬ 
move many instrumental artifacts from the time series 
while preserving any astrophysical variability. We also 
made use of the information in the file headers to ob¬ 
tain estimates of the combined differential photometric 
precision (CDPP). 

We then examined by eye all 2600 O — C curves for 
interesting variations. Some of the results have been 
reported elsewhere (Rappaport et al. 2013; Tran et 
al. 2013). The O — C curve for KIC 9164561 is discussed 
in Sect.J^ We also searched for folded light curves that 
(i) exhibit a flat-bottomed eclipse that is deeper than 
the other eclipse, and (ii) are associated with hot Kepler 
targets, i.e., targets that have a reported Tes > 7000 
K. The deep flat-bottomed eclipse - high temperature 
combination is a signature of a hot white dwarf in or¬ 
bit with a normal, F, A, or B star. In addition to the 
four known Kepler low-mass bloated hot white dwarfs, 
we found two others, previously unreported, in the bina¬ 
ries KIC 9164561 and KIC 10727668. 

The nondegenerate components of KIC 9164561 and 
KIC 10727668 are, according to the Kepler Input Cat¬ 
alog, hot with effective temperatures of ^8000 K and 
10,000 K, respectively, and thus are not susceptible to 
the starspot modulations that are exhibited by cooler 
stars. An examination of the Fourier transform of the 
KIC 10727668 flux time series shows no evidence for 
any periodicities, including those from stellar pulsations, 
other than that of the orbit. In the case of KIC 9164561, 
the Fourier transform of the flux time series shows a set 
of pulsations at frequencies > 20 cycles/day, which does 
not significantly affect our light curve or timing analysis. 
We therefore made no corrections to the data except for 
normalizing each quarter to its median value. 

The orbital periods given in the Kepler binary catalog 
(Matijevic et al. 2012; Kirk et al. 2015) were used to 
produce initial folds of the data. From these folds, we 
established a trial template for the orbital light curve. 
This profile was then used to fit each individual eclipse in 
the time series, thereby determining a set of 893 and 108 
eclipse times for KIC 9164561 and 10727668, respectively. 
We then fit a linear function to the sequence of eclipse 
times of the form: = to + nPorh-, where n is the nth 

eclipse. The two orbital periods thus determined are in 
good agreement with those given in the Kepler binary 
catalog, but have somewhat smaller uncertainties. The 
periods and their uncertainties for these two systems are 
listed in Tabled 

We also fit the eclipse times for each system including a 
quadratic term to allow for possible changes in the orbital 
period or light travel time effects that might occur if the 
binary is part of a triple system. These are discussed in 
Sect. 0 below. 


® https://archive.stsci.edu/kepler/ 
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Fig. 1.— Folded light curve for KIC 9164561 for an orbital period 
of 1.2670400 days. The folded data have been averaged into 3- 
minu te b ins. The red curve is a fit to a model that is described in 
Sect. [3T| The bottom panel shows the residuals of the data minus 
the fitted model. 
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Fig. 2.— Folded light curve for KIC 10727668 for an orbital 
period of 2.305897 days. The folded data have been averaged into 
3-mi nute bins. The red curve is a fit to a model that is described in 
Sect. [3T] The bottom panel shows the residuals of the data minus 
the fitted model. 


3. ORBITAL LIGHT CURVES 

The data for both KIC 9164561 and KIC 10727668 
were folded at the orbital periods we determined from 
fitting the individual eclipses, as described above. The 
results were averaged into 3-minute bins, and are shown 
in Figures [2 and 

Due to the small curvature observed in the eclipse tim¬ 
ing residuals for KIC 9164561 (see Sect.J^ we removed a 
quadratic term from the eclipse phases before producing 
the folded light curve in Fig. The phase corrections 
amount to less than 0.002 days (or 0.0016 orbital cycles). 

In each case, the light curves are normalized to unity 
at the bottom of the deeper eclipse when the hot white 
dwarf has been totally eclipsed. This is then the light 
level of the primary A star at a minimum of its ellip¬ 
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Fig. 3.— Stellar evolution tracks in the mean stellar density, 
(p), and Tejf plane. The stellar masses range from 1.2 to 2.6 Mq 
in steps of 0.2 Mq. The tracks were produced with Yonsei-Yale 
stellar isochrones (Yi et al. 2001). The locations of KIC 9164561 
and KIC 10727668 are indicated by yellow rectangles (1 a confi¬ 
dence regions). The mean stellar densities are determined from the 
eclipse parameters. 


soidal variations. The relatively small separation of the 
two components of KIC 9164561 is manifest by the short 
orbital period, the large amplitude of the ellipsoidal light 
variations (“ELV”), and the duration of the eclipse rela¬ 
tive to the orbital period. 

3.1. Simple Geometrical Model for the Light Curves 

The light curves of these two binaries were analyzed 
with a relatively simple geometrical model which involves 
13 free parameters. In brief, there are 7 free parameters 
associated with the eclipses: i?wd/7?i, Ri/a, i, 
two limb-darkening coefficients, and where subscript 
“1” refers to the primary A star, a is the orbital sepa¬ 
ration, i the orbital inclination angle, and to, the time 
of the middle of a primary eclipse. The out-of-eclipse 
phases are fitted with a set of 6 terms, i.e., a sine and a 
cosine term for each of the orbital frequency and its next 
two higher harmonics. Each of these can be associated 
with different contributions to the ELVs, the illumina¬ 
tion term, and the Doppler boosting effect, as described 
in some detail in Carter et al. (2011). 

Once the best-fitting parameters were found with a 
Levenberg-Marquardt fitting algorithm, the allowed re¬ 
gion in parameter space was defined with a Markov Chain 
Monte Carlo (MCMC) routine, which used the Gibbs 
sampler and the Metropolis-Hastings algorithm, and a 
likelihood proportional to exp(—x^/2) (see, e.g., Holman 
et al. 2006; Sanchis-Ojeda et al. 2013). 

The results for the fitted model are plotted as red 
curves in Eigs. and The high quality of the fit in 
each case is indicated by the white-noise appearance of 
the residuals. The results for the fitted parameters and 
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their uncertainties are summarized in Tabled 


4. THE A STAR COMPANIONS 


The mass and radius of the primary star in each system 
are essential to quantitative models of the two binaries. 
There are two accessible parameters that can be utilized 
to estimate the primary mass: the effective temperature, 
Ti, and the mean stellar density, pi. The Kepler In¬ 
put Catalog (KIC) gives Ti = 8059 K for the primary 
star of KIC 9164561 and 9696 K for the primary of KIC 
10727668. These temperatures are based on 5-color pho¬ 
tometry. For KIC 9164561 we measured = 7870±120 
K from the same spe ctra used for the radial velocity mea¬ 
surements (see Sect. |5.1|). In the case of KIC 10727668, 
where we have no additional spectral information, we 
estimate the uncertainty in the KIC value of Ti at ap¬ 
proximately ± 300 K (Huber et al. 2104). 

The second readily, and more accurately, determined 
parameter associated with each primary star is its mean 
density, pi = Mi/( 47 ri?f/3). This mean density is a 
function only of the light curve-determined parameters 
P and a/i?i, with a weak dependence on q (the mass ra¬ 
tio M 2 /Ml). The dependence of pi on these parameters 
may be found by dividing Kepler’s third law by to 
obtain: 




( 3 ) 


This relation was pointed out by Seager & Mallen- 
Ornelas (2003) in the context of transiting exoplanets. 
For exoplanets, the value of q is much smaller than the 
error in a/Ri so that pi can be determined indepen¬ 
dently of any mass information. While the same is not 
quite true with our two binaries, the value of q is still suf¬ 
ficiently small, i.e., <7 0.1 that a rather accurate value 

of Pi may be determined nearly independently of q. In 
particular, by using Eqn. (3) with the results from our 
photometric analysis of Ri/a and with, as a first approx¬ 
imation, q = 0.1, we find that pi « 0.17 ± 0.04 g cm“^ 
for KIC 9164561, and pi ~ 0.63 ± 0.04 g cm“^ for KIC 
10727668. These values are refined below using estimates 
for q that are based on estimates of the masses of the two 
component stars. 

Figure shows stellar evolution tracks for stars of 
a range of initial masses in the p-TgH plane. The lo¬ 
cations of the two target primary stars are indicated 
by the rectangular boxes. We conclude from this that 
the mass of the primary star in KIC 9164561 must be 
Ml = 2.02 ± 0.06 M 0 (KIC 9164561) and that of the 
primary in KIC 10727668 must be 2.22 ± 0.10 Mq. 


5. RV MEASUREMENT OF KIC 9164561 


The spectroscopic follow-up of KIC9164561 was car- 
ried out with the Fibre-fed Echelle Spectrogra ph (EIES; 
Frandsen & Lindbergl 19991 |Telting et al. 120141) mounted 
at the 2.56-m INordic Optical Telescope (WUT'j of Roque 
de los Muchachos Observatory (La Palma, Spain). Spec¬ 
tra were taken under clear and stable sky conditions at 
11 epochs between May and September 2014. We used 
the 1.3" med-res fibre, which provides a resolving power 


of about 47000 in the spectral range 3600-7400 A. The 
observations taken at each epoch were split into two or 
three consecutive sub-exposures of 1200 seconds each to 


remove cosmic r ay hits. Following the observing strat¬ 
egy described in Buchhave et al. (2010), we monitored 
the radial velocity drift of the instrument by acquiring 
long-exposed (Texp=15-30 sec) ThAr spectra in “sand¬ 
wich mode”, i.e., right before and after each epoch ob¬ 
servation. The data were reduced using a customized 
IDL software suite, which includes bias subtraction, flat 
fielding, order tracing and extraction, and wavelength 
calibration. Radial velocities (RVs) were determined by 
cross-correlating the target spectra with the spectrum of 
Vega, for which we adopted a heliocentric r adial velocity 
of —13.5 kms“^(|Allende Prieto et al.||2004|. 

Table [^lists the FlEb RV measurernents, the signal-to- 
noise (S/N) ratio per pixel at 4650 A, and the full-width 
at half maximum (FWHM) of the FIES cross-correlation 
function (CCF). Full-widths at half maximum were ex¬ 
tracted by fitting a Gaussian to the CCFs between the 
observed spectra and the ATLAS9 model spectrum hav¬ 
ing the s ame photospheric parameters as the target star 
(see Sect. 5.1) and smoothed to the FIES spectral res¬ 
olution. Figure 1^ shows the phase-folded RV curve of 
KIC9164561, based on the phase established from the 
Kepler photometry, and the Keplerian circular-orbit fit 
to the results. The K velocity of the primary A star is 
21.64 ± 0.35 km s“^, and the corresponding 7 -velocity 
is -28.43 ± 0.25 km s'^ 

The corresponding mass function is 


/(Ma) = 0.00133 ± 0.00006 Mq 

For an orbital inclinational angle of 71.6°±0.2° (see Table 
[^, this implies: 

= 0.00155 ± 0.00008 Mq 
(Ml + M2)2 ° 

The mass for the A star of Mi = 2.02±0.06 Mq obtained 
above may be combined with the mass function to obtain 
a white dwarf mass of M 2 = 0.197 ± 0.005Mq. The 
mass ratio for the KIC 9164561 system is therefore q = 
0.097 ±0.004. 

Because KIC 10727668 is faint, i.e., 2.6 Kepler mag¬ 
nitudes fainter than KIC 9164561, we were unable to 
obtain any useful RV measurements for that system. 


5.1. Spectral analysis of KIC 9164561 

We started the spectral analysis by coadding all the 
FIES spectra, after shifting them to a zero-velocity ref¬ 
erence frame, which yielded a S/N of « 55 per pixel at 
4650 A. The spectral parameters of the primary compo¬ 
nent in KIC 9164561 were then derived by fitting this 
spectrum to a grid of synthetic spectra calcu lated with 
the SPECTRU M code (iGray fc Corbally |1994|) using AT- 
LAS9 models (jCastelli V Kurucz||2004 ). We assumed a 
secondary-to-primary tlux ratio of 0.0285, as estimated 
from the eclipse light curve modeling (Table J^. We 
adopted the cali bration equat i ons fr om |Allendb Prieto| 
et al.| (|2004|) and |Bruntt et al.| (|2010|) for the purpose of 
estimating the microturbulent Umicro and macroturbulent 
I'macro velocities. These equations were extrapolated to 
the effective temperature of the star. The spectra l type 
was deterrnined f ollowing the method described in |Gan-| 
dolfi et al.| (2008|). The estimated stellar parameters are 
listed in 'I'able Figure shows the co-added FIES 
spectrum in the spectral region encompassing the Hq,, 
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TABLE 1 

FIES RADIAL VELOCITY MEASUREMENTS OF KIC 9164561 


BJDutc 
{- 2 450 000) 

RV 

kms“^ 

O'RV 

kms“^ 

CCF“ FWHM 
km s“^ 

f^FWHM 

kms“^ 

S/N/pixeF 
@4650 A 

Date of Obser.'^ 
(at start) 

6783.67056 

-9.397 

0.673 

103.6 

1.4 

23 

06 May 2014 

6833.58987 

-38.015 

0.714 

106.7 

1.2 

19 

25 Jun 2014 

6833.67106 

-43.987 

0.768 

102.9 

1.2 

17 

25 Jun 2014 

6844.66674 

-7.539 

0.968 

112.1 

1.2 

11 

06 Jul 2014 

6845.43519 

-38.741 

0.794 

109.1 

1.0 

15 

06 Jul 2014 

6845.71215 

-11.787 

0.690 

102.8 

1.2 

23 

07 Jul 2014 

6901.37427 

-22.004 

0.847 

107.6 

1.2 

15 

31 Aug 2014 

6902.40137 

-43.375 

0.743 

104.6 

1.3 

17 

01 Sep 2014 

6915.38485 

-12.215 

0.784 

106.5 

1.3 

15 

14 Sep 2014 

6915.55901 

-7.337 

0.820 

107.2 

1.1 

14 

15 Sep 2014 

6920.43515 

-14.633 

0.791 

105.6 

1.3 

13 

19 Sep 2014 


(a) Width of the cross correlation function, (b) S/N ratio per pixel at 4650 A. (c) Calendar date in UT at the beginning of the target 
observation. 



Orbital Phase 

Fig. 4.— Radial velocity measurements (points with error bars) 
made with the FIES-NOT spectrometer during the period 2014 
June 15 through 2014 September 15. The uncertainty per point is 
estimated to be ^ 700 m s“^. The curve is the best-fit circular orbit 
with orbital period and phase taken from the Kepler photometry. 
The K-velocity of the A star is 21.64 it 0.35 km and the 7 
velocity is —28.43 km”'^. 


TABLE 2 

Spectroscopic parameters of the primary component of 
KIC9164561 


Parameter 

Value 

Effective Temperature Teff (K) 

7870 ± 120 

Spectroscopic surface gravity log 5 (logio cms“^) 

4.0± 0.1 

Metallicity [M/H] (dex) 

O.Oib 0.1 

Microturbulent velocity t^micro (kms“^) 

2.6 ±0.3 

Macroturbulent velocity t;macro (kms“^) 

11.3 ±0.7 

Projected rotational velocity usinz (kms~^) 

75 ± 5 

Spectral Type 

A7V 


H^, and Mgi b lines, along with the best fitting syn¬ 
thetic model. 


5.2. Constraints on a Third Star in KIC 9164561 

We looked for the spectroscopic signature of a third 
(bound) component to KIC 9164561 by cross-correlating 
the FIES spectra with a set of A- through M-type ob¬ 
served templates. The templates were gathered between 
2011 and 2014, as part of the on-going RV follow-up of 
CoRoT and Kepler stars carried out with FIES@NOT 
(e.g., Gandolfi et al. 2013). 

Assuming that the third star is on or near the ZAMS, 
its mass cannot exceed ^2 Mq - corresponding to an 
A7V spectral type - without significantly violating a 
number of observational constraints (see below). The 
eclipse timing variations detected in the Kepler light 
curve (see Sect.[^ constrain the orbital period of a puta¬ 
tive third star in the system to lie between 8 and 19 years, 
implying that its RV is within ~12 km s“^ of the 7 veloc¬ 
ity of the binary system. The FIES RV follow-up covers 
a time span of about 4.5 months, which corresponds to 
less than 4.6 % of the orbital period of the third com¬ 
ponent, constraining its RV variation to be less than ~4 
km s“^ within the observing window. Because of the rel¬ 
atively high pr^ected rotational velocity of the primary 
star (see Table 0, the almost “stationary” CCF peak of 
the third source is expected to be blended with the broad 
peak of the A7V component. 

The FIES CCFs show neither a peak with constant 
RV superimposed on top of the broad peak of the A7V 
star, nor a skewed profile that would suggest the pres¬ 
ence of a another rapidly rotating staip^ In order to as¬ 
sess the significance of the non-detection and constrain 
the parameters of a possible third body, we built up a 
set of synthetic double-lined spectroscopic binary (‘SB2’) 
spectra using ATLAS9 model atmospheres. We adopted 
Gray’s (2005) calibration for dwarf stars to convert spec¬ 
tral types to stellar parameters. We combined the ro- 
tationally broadened model spectrum of an A7V star - 
resembling the primary component - with that of dwarf 
stars with different spectral types (from A7V to M5V) 
and values of usinf. We added random noise - to ac¬ 
count for the average S/N ratio of the FIES spectra - 

The white dwarf is undetected in our spectroscopic data-set, 
owing to the small secondary-to-primary flux ratio of 0.0285 (see 
Table . 
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Fig. 5.— FIES co-added spectrum of KIC 9164561 (black line) 
encompassing the (upper panel), (middle panel), and Mgl 
b (lower panel) lines. The best fitting ATLAS9 model spectrum is 
overplotted with a thick red line. 


and the contaminant light of the second source - to ac¬ 
count for the white dwarf. We assumed a zero veloc¬ 
ity for the third light and radial velocities in the range 
—21.64 < RV < 21.64 kms“^ for the A7V primary star. 
Finally, we cross correlated the combined SB2 spectra 
with narrow-lined synthetic spectra covering different 
spectra types, in the same fashion as for the observed 
spectra. 

Since the contrast of a CCF peak decreases with in¬ 
creasing usini, we considered the following two cases. 

Spectral type later than F5 V. We used gyro-chronology 
to constrain the usini of the third s tar. Given the age 
of the system (~925 Myr; see Sect. 10.2), the rotation 
period of the third component is expected to be longer 


than about 4 days for spectral later than F5V (Barnes 
et al. 2007). This places an upper limit of 18 km s“^ 
on usinb Although this value likely overestimates the 
true V sin i for the late-type stars, we conservatively as¬ 
sumed 18 km s“^ to be the projected rotational velocity 
associated with the third light. We found that any stars 
with spectral type between K2 V and F5V would be de¬ 
tectable in the FIES CCF (for an example see Fig. |^. 
For later spectral types the peak of the third light is em¬ 
bedded in the noise of the CCF, owing to the relatively 
low S/N ratio of the spectra and the low flux of the third 
star. 



Fig. 6 .— Simulated cross-correlation function (red curve) be- 
tween a narrow-lined KOV synthetic template and the spectrum 
of an hypothetical double-lined spectroscopic binary containing an 
A7V with usini = 75km/s and a KOV star with usini = 18km/s 
(see text for details). The dashed red and blue lines mark the RV 
of the KOV (0 km/s) and the A7V star (21.6 km/s), respectively. 
The black curve, shown for reference, is the same CCF, but without 
the presence of the KOV star. 

Spectral type between F5 V and AIV. For v sin i < 
60 km s“^ the CCF peak of the third component is 
clearly visible on top of the peak of the primary star, 
as described in the previous case. For vsini > 90 km/s 
the CCF would assume a variable skewed profile, whose 
shape changes with the changing RV of the primary star, 
which we do not observe. Finally, for 60 km/s < usini < 
90 km/s the two CCF peaks are merged together and 
larger actual RV variations of the primary component 
would be needed to account for the measured RVs (see 
below). In the case of comparable amounts of light con¬ 
tributed by the A7V star and a putative third star, we 
can say two things. First, the FWHM of the CCF would 
vary in a correlated way with the measured RVs. The 
extra widening can be estimated analytically, and the 
broadening factor is given by [1 + 1/2((5u/(t)^], where 
Sv is the observed (i.e., ‘diluted’) instantaneous value of 
the A7V star RV, and a is the Gaussian width of the 
CCF curve. Since the CCF FWHM remains constant 
at 106.5 km s“^ to within an rms scatter of ^3 km s“^ 
(see Table [^, and a 42 km s“^ shift (i.e., dv = 21 km 
s“^) is expected to widen the CCF by 12 km s“^, which 
is distinctly not observed, we conclude that there is no 
other A- to early F-type star of comparable brightness 
in the system. Second, as mentioned above, a stationary 
line would ‘dilute’ the actual Doppler shift by a factor of 
'^Fa 7 v/{J'A rv + J^s), where Fa 7 v and F 3 are the fluxes 
of the A7V star and the third body, respectively. In this 
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case, the properly inferred mass of the white dwarf (based 
on the RV amplitude) would be a factor of approximately 
(1 + higher. This, in turn, would yield an im¬ 

plausibly high mass for the white dwarf when J-^ exceeds 
•^A 7 V by even a factor of 2. Therefore, the mass of an 
early-type third body is independently constrained to be 
< 2.5 M©. 

Thus, based on our simulations, we can exclude nearly 
any plausible near-MS third star that has a spectral type 
earlier than K2 V. While there is still some small ‘wiggle 
room’ left in parameter space, these constraints likely 
imply an upper limit for the mass of any third star to be 
about 0.8 Mq. 


6 . DOPPLER BOOSTING EFFECT 

In principle, the Doppler boosting amplitude (Loeb 
& Gaudi 2003; van Kerkwijk et al. 2010) can serve as 
a surrogate in place of radial velocity (‘RV’) measure¬ 
ments for measuring the binary K velocity. As formu¬ 
lated by Carter et al. (2011) the Doppler boosting (here¬ 
after ‘DB’) amplitude in a binary is given by: 


Adb = 


Kiai 

c 


a2 F 2 Ml 
oil Fi M2 


(4) 


where the F’s are the fluxes of the two stars in the Ke¬ 
pler band. The a’s are Doppler boosting pre-factors cal¬ 
culated from: 


a = 3 — 


din F^\ 
dlnv / 


(5) 


where F^ is the spectral flux density at frequency 
and a must be evaluated using stellar atmosphere mod¬ 
els averaged over the Kepler band. The second term in 
Eqn. Q represents the contribution to the DB ampli¬ 
tude from the white dwarf in the system (note the minus 
sign). This method for evaluating the K velocities has 
been verified by Ehrenreich et al. (2011) and Bloemen et 
al. (2012) for the case of KOI 74, a system quite analo¬ 
gous to the ones discussed in this paper. 

KIC 9164561: The measured Doppler boosting am¬ 
plitude in this system, i.e., the amplitude of the sinojt 
term, where uj is the orbital frequency, and t is time, is 
-10 ± 26 parts per million (‘ppm’). Given the Ki ve¬ 
locity obtained from the RV measurements (see Sect.[^, 
for nominal values of ai, Adb for this system should oe 
«160 ppm. One possible explanation for this discrep¬ 
ancy is based on the large illumination semi-amplitude, 
i.e., the coefficient of the cos wt term, that is measured to 
be 2385 ± 36 ppm. In this case, if the phase of the illu¬ 
mination is shifted in the proper sense from its nominal 
alignment with the eclipses by as little as 4°, the en¬ 
tire DB effect would be cancelled. We note that in KOI 
74 the illumination term is 600 times smaller because 
of the wider orbital separation and the lower luminos¬ 
ity of the white dwarf, and is therefore very unlikely to 
materially affect the observed Doppler boosting ampli¬ 
tude. In order for the phase of the illumination effect 
in KIC 9164561 to be shifted by the amount required to 
cancel most of the Doppler boosting signal, the A-star 
primary would have to be rotating faster than the rate 
that is synchronous with the orbit, such that some of the 
heated surface material is carried into the night side of 
the star. The measured rotation velocity of this star is 
« 75 ± 5 km s“^, a value that is to be compared with 


the equatorial rotation velocity of a star of radius 2.5 
i ?0 (see Table rotating with a period of 1.267 days, 
i.e., « 100 km s“^. Thus, if anything, the primary is ro¬ 
tating subsynchronously and the DB signal is therefore 
expected to be enhanced relative to its nominal ampli¬ 
tude. This is a potentially important discrepancy which 
merits follow-up study. 

KIC 10727668: In this system, the Doppler boosting 
signal is measured to have a significantly non-negative 
value, 143 ± 20 km s“^. The illumination amplitude, 
557 ppm, is 4 times smaller than in the case of KIC 
9164561; this system may therefore be less susceptible to 
an illumination-induced perturbation of the DB ampli¬ 
tude. If Eqn. is used to evaluate Ki we find: 

20.0 ±2.8 , _i 

= (1-0.02M./M,) * <'’> 

where we have taken ai = 2.15, 012 = 1-78 (van Kerkwbk, 
private communication), and F 2 /F 1 ~ 0.025 (see Fig. 1^. 
This expression for Ki may be entered into the mass 
function, and the cube root of the result taken to obtain 
the following equation involving the masses of the two 
stars: 

M 2 _ 0.107 ±0.015 1/3 

{Mi+M2f^ (1 - O.O 2 M 1 /M 2 ) ® 

where we have taken the orbital inclination angle to be 
86.4° (see Table |^. This equation can be solved nu¬ 
merically for M 2 given any specific value for Mi. For a 
primary star mass of 2.2 M© (see Fig. [^, the mass of 
the secondary is 

M 2 ~ 0.266 ± 0.035 M© 

where the error bar takes into account the uncertainty in 
Ml as well as the uncertainty in Ki. Again we caution 
that this estimate for M 2 does not take into account any 
systematic error due to an “amplitude leak” in the DB 
term. 


7. ELLIPSOIDAL LIGHT VARIATIONS 


The ellipsoidal light variations (‘ELVs’) seen in the out- 
of-eclipse region of the light curves also contain informa¬ 
tion on the mass ratios, q = M 2 /Mi. The basic expres¬ 
sion for the ELVs (Kopal 1959; Morris & Naftilan 1993; 
Carter et al. 2010) is: 


and 


Belv - -Zi{2) q {Ri/aY sin^ i 


Zi{2) 


45±3u 
20(3 - u) 


(1 ±t) 


( 8 ) 

(9) 


where i?ELV is the amplitude of the cosine term at twice 
the orbital frequency, u is the linear limb-darkening coef¬ 
ficient, and r is the gravity darkening coefficient. For the 
hot primary stars in KIC 9164561 and KIC 10727668 we 
take u = 0.3 and r = 1. This leads to a numerical value 
for Zi{2) of 1.68. We can then utilize the measured val¬ 
ues of Belv of —8880 ppm and —1179 ppm for the two 
systems, respectively (see Table Similarly, we take 
the values for Ri/a and i for eacn system from Table 
Finally, we solve for q to find: 

gELv(9164561)= 0.104 ±0.002 (10) 

<7elv(10727668)= 0.15 ±0.01 (11) 
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where the error bars are statistical only and do not in¬ 
clude any systematic uncertainties in the use of Eqn. (|^ 
to evaluate q from the ELV amplitude. 

The value of q found for KIC 9164561 from the RVs 
is 0.097 ± 0.004 (Sect. 15.11) which is in quite satisfactory 
agreement with ^elv- 

On the other hand, the value of q for KIC 10727668 
found from the ELV analysis (Eqn. ([Tl])), while still in 
marginal agreement with the value obtained from the 
Doppler boosting analysis (0.12 ± 0.02), would imply a 
white dwarf mass as high as M^d — 0.33 Mq. This is 
quite on tlw high side of what we expect for this system 
(see Sect. 10). 


8 . MASSES AND RADII OF THE WHITE DWARFS 


In order to determine the sizes of the binary systems 
we need the total masses. In the case of KIC 9164561, 
where a direct RV measurement has been made, the to¬ 
tal system mass is 2.22 ± 0.07 Mq. With the appli¬ 
cation of Kepler’s third law, this leads immediately to 
an orbital separation of a = 6.40 ± 0.077?©. For KIC 
10727668 the Doppler boosting analysis yielded a mass 
ratio q ~ 0.12. If we tentatively accept this mass ratio, 
and apply Kepler’s third law, we find the orbital sepa¬ 
ration is a = 9.95 ± 0.15i?©. The reason that the un¬ 
certainty is quite small, in spite of the larger uncertainty 
in q for this system, is that a depends only on the cube 
root of the total system mass, and g is a small number. 
These values for the orbital separation are summarized 
in Table [3] 

The radii of the primary and secondary stars follow im¬ 
mediately from the Ri/a and i? 2 /a values derived from 
the light curve analyses and the values of the orbital 
separations. These radii are summarized in Table 
Of particular note are the large radii of the two white 
dwarfs, 0.277 i?© and 0.151 i?© for KIC 9164561 and 
KIC 10727668, respectively, which indicate that these 
stars are quite thermally bloated. The degenerate radii 
for cold white dwarfs of either He or hybrid HeCO com¬ 
positions, are given by Eqn. (4) in Nelson & Rappaport 
(2003). The degenerate radii for the white dwarfs in KIC 
9164561 and KIC 10727668 are 0.020 i?© and 0.018 i?©, 
respectively, implying thermal bloating factors of 14 and 
8 . 


9. EVIDENCE FOR A THIRD STAR IN KIC 9164561 

We carried out an eclipse timing analysis for both bina¬ 
ries discussed in this paper. We used a code to produce 
the observed minus calculated curve {‘O — C’ curve) that 
was tested in our search for triple star systems (Rappa¬ 
port et al. 2013). Only the O — C curve for KIC 9164561 
showed an interesting non-linear behavior. The results 
are shown in Fig. with black points for the average of 
the O — C curves for the primary and secondary eclipses. 
The O — C curve shows two notable features: (i) a clear 
curvature, and (ii) small-amplitude oscillations with a 
period of ^170 days. The small amplitude oscillations 
are due to an artifact that results from a beat between 
the orbital period and Kepler long-cadence sampling (see 
Eqns. (4) and (5) in Rappaport et al. 2013). 

The quadratic behavior in the O—C curve is highly sig¬ 
nificant (35 a), and can be interpreted as either a change 
in the orbital period, or as a Roemer delay due to the 
motion of the binary induced by a third body. The Hrst 



Fig. 7.— O — C curve averaged over the primary and secondary 
eclipses of KIC 9164561 plotted as black points. The O — C curve 
is based on a binary period of 1.267040 days (see Table . It is 
evident that the O — C curve has a distinct curvature - significant 
at the 35-0- statistical level. In addition, there is a small-amplitude 
oscillation with a predictable period of ~170 days which is an ar¬ 
tifact caused by a beat between the orbital period and the Kepler 
long-cadence sampling time. The characteristic timescale for the 
orbital period change implied by the overall long-term quadratic 
term is ~2 x 10^ years. If, instead, the curvature is attributed 
to a Roemer delay induced by a third body, the inferred orbital 
period is in the range of 8-19 years. Here the red curve is the best¬ 
fitting circular orbit with a period of 13.7 years, and it includes 
a sinsuoidal modulation with a 170-day period to account for the 
artifact caused by the beats. 


option requires changes in the binary orbital period on 
a characteristic timescale of P/P ~ 2 x 10^ years. Since 
there is no significant mass transfer taking place in the 
system at the current time, this cannot be the cause of 
the apparent orbital period changes. One would have to 
invoke another mechanism, such as the Applegate effect 
(Applegate & Patterson 1987; Applegate 1992; Kotze & 
Charles 2010) to explain the short timescale changes in 
the orbit. We focus here on a more interesting and, in 
our view, more likely scenario, i.e., that there is a third 
body in the system. It is then possible to constrain the 
orbital radius of such a third body as well as its mass, if 
we restrict our attention to coplanar and circular orbits. 
Otherwise, the possibilities are nearly limitless. 

The results of such constraints are shown in Fig.[^ We 
find that the orbital period must lie in the range of 8 < 
Ttrip < 19 years in order to provide at least a somewhat 
satisfactory fit to the O — C curve (±4 ct limits). The 
blue curve shows the best fit (projected) orbital radius 
as a function of Ptrip- The red curve is the corresponding 
mass function, while the green curve is the mass of the 
third body under the assumptions that the mass of the 
inner binary is 2.2 M© and that the orbital inclination 
of the third body orbit is 90°. Since the mass of the 
third body, if it is assumed to be on or near the ZAMS, 
cannot exceed ~2 M© without affecting the spectrum, 
RV analysis, and the third ligh t of the syst em in ways 
that are not observed (see Sects. 5.1 and 5.2), the orbital 
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Fig. 8. — Illustrative constraints on the projected orbital radius 
(blue curve), mass function (red curve), and mass of the third body 
(green curve) in KIC 9614561 based on the O — C curve shown in 
Fig. Only third-body orbital periods between 8 and 19 years 
yielcisomewhat acceptable fits (marked by dashed vertical lines), 
and periods greater than 14 years are excluded because the inferred 
third-body would be too massive and therefore likely too luminous 
to be missed in our spectroscopy. Only circular coplanar orbits 
were considered in setting these constraints. 

period is further restricted to the range of 8 < Ptrip < 14 
years. If we adopt the formal constraint on the mass of 
a main-sequence third star in the system from Sec. |5.2| 
of M 3 < 0.8 Mq we find that: 0.5 < M 3 < 0.8 Mq and 
8 ^ ^3 ^ 10 years. 

To this point we have not yet considered the possibility 
that the putative third body might be a white dwarf. The 
orbital period would then be restricted to the range 8—12 
years since the mass would have to be < 1.4 Mq (see 
Fig.§. If there had been a prior episode of mass transfer 
from the progenitor of the outer body white dwarf to the 
inner binary (KIC 9164561; see, e.g., Tauris & van den 
Heuvel 2014) then, according to Eqn. Q, the mass would 
lie in the range of 0.6 — 0.8 Mq . Such -^ite dwarfs would 
have cooled to the point where they would contribute 
very little light to the system (< lO^'^Ltot)- 

Based on the inferred luminosity of KIC 9164561, and 
the appropriate bolometric corrections, we infer that its 
distance from the Earth is « 390 pc. Typical orbital 
separations to the putative third body, for periods in the 
range of 8—19 years, are expected to be « 5—10 AU. At a 
distance of 390 pc, the corresponding angular separations 
are then 14—26 mas. This may be detectable with Gaia. 

10. DISCUSSION 

Here we discuss the properties of the KIC 9164561 
and KIC 10727668 systems and some of the evolution¬ 
ary paths that may have led to their formation. We 
can then check whether the expected system parameters 
match those determined empirically. 

10 . 1 . Evolutionary Pathways 


TABLE 3 

KIC 9164561 AND KIC 10727668 System Parameters 


parameter 

KIC 9164561 

KIC 10727668 

RA (J2000) 

19h 42m 27.64s 

19h 21m 45.58s 

Dec (J2000) 

45° 30' 17.0" 

48° 02' 42.4" 

P [days] 

1.2670400 ±0.0000002 

2.305897 ± 0.000007 

to [MJD]'^ 

64.61956 ±0.00013 

1207.3677 ±0.0005 

Kp [mag.] 

13.71 

16.32 

R 2 /R 1 

0.1091 ±0.0016 

0.0893 ± 0.0006 

Rl/a 

0.397 ±0.003 

0.170 ±0.003 

a [Pq] 

6.40 ±0.07 

9.95 ± 0.15 

Rk,2IRk,G 

0.0285 ±0.001 

0.0245 ±0.001 

inclin. [deg] 

71.59 ±0.22 

86.45 ±0.5 

Ml [Mq] 

2.02 ± 0.06 

2.22 ± 0.10 

Rl [.R©] 

2.54 ± 0.03 

1.69 ± 0.04 

Ti [K] 

7870 ± 150 

9696 ± 300 

Li [Lq\ 

22.3 ± 1.5 

22.9 ± 3.1 

Pl [g/cc] 

0.171 ±0.04 

0.63 ± 0.04 

log ( 91 ) [cgs] 

3.93 ± 0.02 

4.32 ± 0.04 

M 2 [Mq] 

0.197 ± 0.005 

0.266 ± 0.035 

R 2 [P©] 

0.277 ± 0.005 

0.151 ± 0.004 

T2 [K] 

10410 ± 200 

14110 ± 440 

L 2 [I/©] 

0.82 ± 0.07 

0.82 ± 0.11 

P 2 [g/cc] 

12.8 ± 1.2 

120 ±19 

log ( 92 ) [cgs] 

4.84 ± 0.04 

5.55 ± 0.06 

Harmonic Amplitudes [ppm] 

Al [sin((/))]'= 

-10 ± 27 

143 ± 20 

Pl [cos(()4)]‘^ 

2385 ± 36 

557 ± 22 

A 2 [sin(2()i)] 

-602 ± 27 

-47 ± 20 

P 2 [cos(2</))]® 

-8880 ± 35 

-1179 ± 21 

A 3 [sin(3?i)] 

-30 ± 28 

-4 ± 19 

P 3 [cos(3</))] 

-888 ± 36 

-48 ± 23 


Note. — a = BJD-2454900; b = ratio of white dwarf flux to 
primary A-star flux in the Kepler band; c = Doppler boosting 
amplitude; d = main illumination term; e = main ellipsoidal light 
variation term. See §|^and §[3 for details. 

There are two important evolutionary pathways that 
could lead to the formation of low-mass helium white 
dwarfs in close binaries (see, e.g., van Kerkwijk et al. 
2010; Carter et al. 2011; Breton et al. 2012). The first 
channel includes one episode of common envelope (CE) 
evolution wherein the more massive star of the primor¬ 
dial binary overflows its Roche lobe after having become 
a giant with a degenerate core, and then undergoes dy¬ 
namically unstable mass transfer onto its binary com¬ 
panion (see, e.g., Taam et al. 1978; Webbink 1984; Taam 
and Bodenheimer 1992). The rapid rate of mass transfer 
causes the companion to spiral inward towards the de¬ 
generate core of the giant while concomitantly expelling 
the envelope of the giant into the ISM. The resulting bi¬ 
nary is composed of a degenerate white dwarf in a tight 
circular orbit with a largely unchanged companion star. 

The second channel involves stable mass transfer be¬ 
tween the components of the progenitor binary. Unfor¬ 
tunately the degree to which mass transfer is conserva¬ 
tive and the effectiveness of orbital angular momentum 
loss mechanisms (e.g., magnetic braking; see, for exam¬ 
ple, Verbunt & Zwaan 1981) are uncertain. For reason¬ 
able physical assumptions in regard to the mass transfer 
and angular momentum loss processes, this channel can 
explain the formation and observed properties of these 
two systems (see Nangreaves et al. 2015, and references 
therein). 
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(i) The Common Envelope Channel 
Although a physically realistic 3-D hydrodynamic sim¬ 
ulation of CE evolution has not yet been carried out, 
a reasonably good approximation of the effects of the 
process can presumably be obtained from first principles 
using energy conservation. As the companion star moves 
at high velocity through the envelope of the giant, a frac¬ 
tion ^ of the orbital energy is transferred to the tenuous 
envelope material. The companion star then spirals in¬ 
ward and more and more energy and angular momentum 
are transferred to the giant envelope. The gravitational 
potential energy that is released may be sufficient to eject 
the envelope. It is assumed that this process occurs on a 
timescale much shorter than the thermal timescale of the 
giant star. One of the more commonly used CE energy 
constraints was analyzed by Dewi & Tauris (2000). It is 
given by 


G (Me + Me) Me 


XR 


1,0 


= e 


GMeM2.0 G(Me + Me)M2.0 


2 a/ 


2ao 


( 12 ) 

where Me, and M 2 ,o are the core and envelope 
masses of the WD progenitor, and the mass of primordial 
secondary, respectively (see, e.g., Webbink 1984; Pfahl 
et al. 2003). Note that ^ is a dimensionless measure 
of the efficiency of the conversion of orbital energy that 
goes toward unbinding the envelope of the giant (the en¬ 
ergy required to eject the envelope, AEb, must satisfy 
AEb = ^AAorb) and is usually taken to be < 1. The 
dimensionless factor A is intended to take into account 
the actual binding energy of the giant at the onset of the 
CE phase. It largely depends on the ratio of the core 
mass to the pre-CE-phase mass of the giant. Values of 
A have been tabulated for many different evolutionary 
states and masses of giants by Tauris & Dewi (2001) and 
Podsiadlowski et al. (2003b). For the types of giants that 
could produce helium white dwarfs, A is of order unity, 
normally ha ving a value between 0.5 and l.Op^ Also note 
in equation (12) that Me, the core mass ofthe giant, is 
approximately equal to the final white dwarf mass, M^d- 
Thus the original mass of the primary can be written as 
Ml 0 = Me -I- Me and its radius is given by Rip. Finally, 
a/ is the separation of the white dwarf and secondary 
star in their mutual circular orbit after the CE phase 
while oq is the separation of the components of the pro¬ 
genitor binary. Using Kepler’s third law, equation (12) 
can be rewritten as 


qi (l±^\_ J (^\ ^ ^ 
90 VI+ 9// V^// J V ao y 

where ry = ^A is a dimensionless parameter that is 1. 
Also qo{= Mip/M 2 p) and <//(= Mc/M 2 ,o) are the mass 
ratios of the primordial and present day binaries, respec¬ 
tively. At the onset of (unstable) mass transfer, the ratio, 
Rip/ciQ, is approximately equal to the Roche lobe radius 
of the progenitor primary divided by the primordial sep¬ 
aration; its value is constrained by Roche geometry and 
can be expressed as a simple function of go (see, Eggle- 
ton 1983). Given that the radius of the red giant will 
have a very strong dependence on M^ and a very weak 


V 

2(90 - 9/) 


If we take A < 0.5, which is possible in certain circumstances, 
this only makes the argument presented here more robust. 


dependence on Mip, the constraint of Roche geometry 
allows the approximation of Pq by a function of Me (see, 
e.g., Rappaport et al. 1995). 

Assuming a specific value of rj, and using the observed 
orbital period of the binary (P/) and the inferred mass of 
the present-day companion Mi (= M 2 J = M 2 p), strong 
limits can be imposed on the possible values of the white 
dwarf mass (we take M^d = Me + O.OIMq), by solving 
Eqn. (13) implicitly for M^d and applying the constraint 
that Mip > M 2 P (i.e., go > !)• K is reasonable to sup¬ 
pose that rj will lie in the range 0.5 to 1.0. Taking the 
upper limit of this range for rj, we find that the only 
acceptable values of the core mass for KIC 9164561 are 
Me ~ 0.36 - 0.38 Mq while for KIC 10727668 they are 
Me — 0.40 — 0.42 Mq. Assuming the lower limit of this 
range of ry, we find that the only acceptable values of the 
core mass for KIC 9164561 are Me — 0.39 — 0.40Mq, 
while for KIC 10727668 they are Me ~ 0.44 - 0.45 Mq. 
These ranges of value^^ of Me are clearly incompatible 
with the estimated values of M^d for both KIC 9164561 
and KIC 10727668 based on the cooling evolution of the 
white dwarfs (see the next section), as well as with the 
empirically determined value for KIC 9164561 from RVs. 
We conclude that for a reasonable range of common enve¬ 
lope efficiencies, the formation of either or both of these 
systems via the CE channel is rather unlikely. 

(ii) The Algol Channel 

The other channel involves stable mass transfer during 
the evolution of the progenitor binary. In this scenario, 
the star that is more massive (the donor) evolves more 
quickly and fills its Roche lobe on its nuclear timescale 
before its less massive companion (the accretor) can do 
so. Even if the mass ratio is as large as go — 2, it is 
entirely possible for the mass transfer to proceed sta- 
bl}p^ if the donor has a substantial radiative envelope 
(Ndson & Goliasch 2014). As an example of evolution 
that could explain the formation of the two observed sys¬ 
tems, consider a 1.5 Mq donor that has evolved beyond 
the TAMS and lies in the Hertzsprung gap when it fills 
its Roche lobe. This star can stably transfer mass to a 
1.2 Mq companion on approximately a nuclear and/or 
thermal timescale (Case AB evolution). At this mass 
transfer rate (of ~10“®Mq yr“^), the accretor is not 
likely to expand and hll its own Roche lobq^ During 
the mass transfer process, the donor can evolve up the 
red giant branch and form a substantial degenerate he¬ 
lium core. Since the radius of the giant has a very strong 
dependence on Me but a very weak dependence on its 
total mass, the giant will continuously expand with a 
concomitant increase in the core mass until the envelope 
surrounding the core is so depleted (to a mass of the or¬ 
der of 0.01 Mq; the remainder having been transferred 
to the companion star) that it cannot be sustained and 
ultimately collapses onto the remnant degenerate helium 
core. If approximately 1 Mq is transferred conservatively 
to the accretor during this process, the resulting system 
may then resemble either of the observed ones. Details 


The lower limits arise from the requirement that go > li by 
definition, while the upper limits correspond to go 5 1-2 to take 
into account the fact that we are considering only primaries of 
masses < 2.5 Mq so that only He WDs are produced. 

Mass transfer could then proceed on a thermal timescale. 

If it were to fill its own Roche lobe then this type of evolution 
would probably lead to either a contact binary or a merger. 
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concerning the possible progenitor evolutions will be dis¬ 
cussed in a future paper (Nangreaves et al. 2015). 

10 .2. Cooling Evolution 

Assuming that stable mass transfer from a red giant 
led to the formation of the bloated helium WDs in both 
KIC 9164561 and KIC 10727668, it is possible to check 
whether the observed properties of the WDs are consis¬ 
tent with cooling models of degenerate helium dwarfs. 
Moreover, if the properties of either WD are found to be 
consistent, it is then possible to use it as a chronometer 
in order to determine the age of the binary after the ter¬ 
mination of mass transfer, i.e., just after the collapse of 
the giant to form a helium WD (HeWD). 

The study of the formation and cooling of HeWDs from 
progenitor red giants has been carried out by a number 
of groups (see, e.g., Podsiadlowski et al. 2003; Nelson et 
al. 2004; Althaus et al. 2005, 2013; Istrate et al. 2014). 
We use the code described by Nelson et al. (2004) to ana¬ 
lyze the cooling evolution of very low-mass helium WDs. 
It is assumed that the giant (i.e., the WD progenitor) 
had a solar metallicity {Z = 0.02), and that its initial 
mass was between 1.0 and 1.5 Mq. The initial hydro¬ 
gen abundance (by mass) of the progenitor star while 
on the ZAMS is taken to be equal to 0.71 and the ratio 
of the mixing length to the pressure scale-height ratio 
{£/Hp) is set equal to 1.5. The initial orbital separation 
was carefully selected so that only very low mass helium 
WDs were created. Diffusion (gravitational settling and 
radiative levitation) was not included in the computa¬ 
tions because diffusion is a very fragile physical process 
that can be negated by turbulence in the outer layers of 
these rapidly rotating (proto-)white dwarf stars. More¬ 
over, the inclusion of diffusion has little effect on the 
cooling times for the models of interest (i.e., bloated he¬ 
lium WDs) because these systems are very young and 
have not even evolved to reach their maximum surface 
temperatures before joining the usual WD cooling track 
in the HR diagram. It is only after the WDs have passed 
through their maximum Teff that the compression of the 
outer hydrogen-rich layers, due to contraction, can cause 
the residual hydrogen to start to burn, leading to a rapid 
expansion of the outer layers as a result of a thermonu¬ 
clear runaww (TNR). 

In Figure we show the results of the cooling evolution 
of several models that have masses near or smaller than 
the TNR lower mass limit, i.e.. Me < 0.21 Mq. The evo¬ 
lution of the radius of the WD is shown as a function of 
its effective temperature, T^s- When the envelope of the 
progenitor giant collapses, the radius of the proto-WD 
continually decreases, as long as no TNRs occur, until it 
asymptotically reaches its zero-temperature radius. The 
surface temperature of the object initially increases (this 
increasing temperature defines the proto-WD phase) un¬ 
til it reaches a maximum value of where it then fol¬ 
lows a cooling track and asymptotically approaches zero 
temperature. 

Assuming that the white dwarfs in both KIC 9164561 
and KIC 10727668 formed according to the scenario dis¬ 
cussed above, we can infer reasonably precise values of 
the mass based on the observed radii and values su¬ 
perposed on the cooling tracks (see Fig. [^. We find that 
the mass of the WD in KIC 9164561 is 0T89± 0.010 Mq, 
and that of the WD in KIC 10727668 is 0.201±0.010Mq. 



Fig. 9. — Cooling tracks in the radius— Tgff plane for low-mass 
white dwarfs that have lost their outer envelopes via mass transfer 
to a companion star. The tracks are labeled by the mass of the 
white dwarf. The filled red circles are the approximate current lo¬ 
cations of the six known thermally bloated hot white dwarfs found 
among the Kepler binaries. The open blue circles are eight low- 
mass white dwarfs found in double-degenerate systems (Hermes et 
al. 2014), while the open green circles are for systems similar to 
those reported here, and that were found using the WASP photo¬ 
metric database (Maxted et al. 2013; 2014); these have been added 
for comparison. 


These masses agree with those inferred from the cooling 
evolutions of Althaus et al. (2013) to within 5 to 10% 
(their values of the mass being smaller), and to within 
< 2% for the values inferred from the work of Istrate 
et al. (2014). We also carried out this same type of 
WD cooling calculations with the MESA code (Paxton et 
al. 2011) and found very similar results. The ages of these 
objects as measured from the cessation of mass transfer 
are approximately 925 Myr and 950 Myr, respectively. If 
the companion star in KIC 9164561 has a mass of 2.0 Mq, 
then the time required to reduce its average density to 
the inferred value of 0.17 g cm“^ is in good agreement 
with the theoretical predictions for the evolution of a so¬ 
lar metallicity star. The orbital period predicted by our 
models is approximately 1.7 days, but it should be noted 
that our assumptions about the degree of mass conserva¬ 
tion during the mass-transfer phase and the magnitude of 
angular momentum losses leading to the formation of the 
giant are not precisely known, and this uncertainty could 
have a significant impact on the final value of the orbital 
period. As for KIC 10727668, the predicted orbital pe¬ 
riod from our models almost matches the observed one. 
However, the predicted age and somewhat higher mass 
of the companion star are very problematic because the 
observationally inferred average density implies that the 
companion star has not significantly evolved from its ini¬ 
tial position on the ZAMS. Based on the implied cooling 
age of the system and a mass of 2.2 Mq, we would ex¬ 
pect the star to be much closer to the TAMS with a 
correspondingly lower density (contrary to what is im¬ 
plied in Figure 3). If the mass of the companion were 
substantially less (say, R:!l.7 Mq), then a WD composed 
primarily of helium and with a mass of 0.21 Mq would 
be entirely compatible with the observations. 

Based on the cooling evolution, we conclude that the 
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theoretical models for the evolution of KIC 9164561 lead 
to properties that agree well with the observationally in¬ 
ferred properties. The predictions for KIC 10727668 are 
problematic, especially given the implied cooling age of 
the white dwarf in the system. For an age in the range 
of 100-200 Myr, a white dwarf mass of 0.27 (in¬ 
ferred from the Doppler boosting effect; see Sect.|^, and 
Teff ~ 14, 000 K, it is very difficult to model the proper¬ 
ties with any type of solar metallicity HeWD. To under¬ 
stand why this is true, consider the evolutionarily limit¬ 
ing case for which all of the hydrogen has been stripped 
from the atmosphere and envelope of the HeWD. Based 
on the WD cooling tracks of Nelson & Davis (2000), a 
0.27 Mq WD of this type can reach a temperature of 
14,000K at either of two ages: (i) considerably less than 
10 Myr; or, (ii) « 250 Myr. The luminosity of the WD 
for the latter case is two orders of magnitude smaller 
than what is inferred (see Table and can be auto¬ 
matically ruled out. For the other solution (correspond¬ 
ing to the early-WD phase before the maximum is 
reached), the luminosity is approximately correct but a 
cooling age of only several million years after the birth of 
the WD seems to be unacceptably small based on prob¬ 
ability arguments. If a thick layer of hydrogen were to 
be added to the envelope, the cooling would be slowed 
but incompatibilities with the observations would persist. 
For example, the luminosity of the HeWD while it is still 
young would be considerably higher and it would likely 
undergo a TNR. Perhaps we are observing the HeWD 
while it is undergoing one of its TNR cycles; the prob¬ 
ability of observing it at this epoch may be unaccept¬ 
ably small. Another possible explanation might be a 
difference in the actual metallicity from that which we 
assumed; however, making the metallicity smaller gen¬ 
erally makes the incompatibilities larger. Thus if the 
mass of the WD in KIC 10727668 is actually about 0.27 
Mq , then it might have been derived from a higher mass 
progenitor (> 2.5 Mq) or perhaps via a different evolu¬ 
tionary pathway than we have considered here. We plan 
to investigate all of these possibilities in a future paper. 

For the sake of completeness, we have included the ro¬ 
bust measurements of effective temperatures and radii of 
eight other eclipsing or tidally distorted, extremely-low 
mass WDs (found in double degenerate systems; see Her¬ 
mes et al. 2014). The observations are plotted in Fig. 
as open blue circles. It should be noted that these addi¬ 
tional WDs have properties that are reasonably similar 
to the six hot WDs found by Kepler except that they are 
typically less bloated. Based on the cooling models, we 
expect the masses of all of these systems to reside in the 
range of 0.15-0.25 Mq, which is generally consistent with 
the values reported in Hermes et al. (2014). Finally, in 
this regard, we note that these white dwarfs have a very 
different evolutionary scenario from the systems exam¬ 
ined in our work in that they are likely the result of two 
common-envelope episodes. 

10.3. Future Evolution 

Here we examine very briefly what the future holds 
for systems such as KOI 74, KOI 81, KOI 1224, KOI 
1375, KIC 9164561, and KIC 10727668. When these sys¬ 
tems evolve to the point where the primary A star fills 
its Roche lobe and commences mass transfer onto the 
white dwarf, their He core masses will range from 0.19 


to 0.27 Mq for the shortest (KIC 9164561) and longest 
(KOI 81) orbital periods of these systems, respectively 
(as estimated using Eqn. .<§)■ Once the mass transfer 
onto the white dwarf begins, it will be dynamically un¬ 
stable due to the very large mass ratio (typically 10:1), 
and a common envelope will ensue. Using Eqn. ( |12[ ), we 
can compute an approximate expression for the ratio of 
the final orbital separation to the initial orbital separa¬ 
tion, and we find: 


o/ ^ MwdAle iirLiq) 
a, ~ Me MA 2 


(14) 


where Mwd is the mass of the hot white dwarfs being 
studied in these systems. Mg and Me will be the envelope 
and core mass of the A star when it evolves to fill its 
Roche lobe, Ma = Me + Mg, rj is the common envelope 
parameter of order unity discussed above, and is the 
Roche-lobe radius of the A star in units of the orbital 
separation (with an exp ected numerical value of « 0.6). 

If we apply Eqn. (141 to estimate the final orbital sep¬ 
aration that will result from this upcoming common en¬ 
velope phase, we find that af F 0.003a^. The current 
orbital separations of the six bloated hot white dwarfs 
in the Kepler field range from 6.4 to 48 i?©. From this 
we can infer that the final orbital separations will be 
a/ < 0.02-0.14 R©, for KIC 9164561 to KOI 81, respec¬ 
tively. 

All of these will be binaries comprised of two roughly 
equal mass white dwarfs each of « 0.2 Mq, if the core of 
the A star does not merge with the initial low-mass white 
dwarf in the system. The corresponding orbital periods 
and Roche-lobe radii would range from 3/4 to 14 min, 
and 0.007 — 0.05 7?©, respectively. All but the largest of 
these Roche-lobe radii are too small to harbor the white 
dwarf cores that will be unveiled by the common enve¬ 
lope. Thus, we conclude that the orbital periods of the 
group of double degenerate systems recently studied by 
Brown et al. (2014) and Hermes et al. (2014) result from 
much wider systems. This is emphasized by the fact that 
the masses of the primary in the double degenerate bi¬ 
naries in the latter studies are typically much more mas¬ 
sive than would result from the six binaries containing 
bloated hot white dwarfs in the Kepler field. 


11. SUMMARY AND CONCLUSIONS 


In this study we report two new low-mass, thermally 
bloated, hot white dwarfs in binary systems with com¬ 
panion A stars. The systems are similar in their basic 
nature to those in KOI 74, KOI 81, KOI 1224, and KOI 
1375. This brings to six the total number of such sys¬ 
tems in the Kepler field. Of the six systems, the two new 
systems have the shortest orbital periods at 1.27 d and 
2.3 d, while KIC 9164561 has the largest, i.e., the most 
thermally bloated, of the six white dwarfs. Its white 
dwarf has a radius that is 14 times larger than its fully 
degenerate radius—it is more than 1/2 i?© in diameteii 
The mass of the white dwarf in KIC 9164561 is 0.197 
± 0.005 Mq and is determined from direct RV measure¬ 
ments of the primary A star as well as some inferences 
about the mass of the A star (i.e., the mean density and 
Teff lead to an estimate for Mi). By contrast, the mass 
of the white dwarf in KIC 10727668 is more uncertain at 
0.27 ± 0.03 Mq. Since there are no RV measurements 
of this fainter system, we relied on measurements of the 
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Doppler boosting effect to determine the orbital veloc¬ 
ity of the A star. The accuracy of this latter mass is 
not quite sufficient for detailed comparisons with binary 
evolution or white-dwarf cooling models. We strongly 
advocate that a larger telescope be used to determine an 
accurate RV for this system. 

We have shown that KIC 9164561 is very plausibly 
part of a triple-star system with an outer period of > 8 
years. The mass of the third star is likely constrained to 
be < 0.8 M©. 

Binary evolution scenarios involving two low-mass, 
e.g., « 1.1 and 1.5 M©, primordial stars that undergo 
stable, conservative mass transfer can nicely explain the 
current mass and evolutionary state of both the compan¬ 
ion star and the hot white dwarf in KIC 9164561. For 
this system, we estimate an age of ~925 Myr since the 
envelope of the white dwarf progenitor was transferred 
to the current A star. By contrast, it is more difficult 
to explain the current state of both the A star and its 
white dwarf companion in KIC 10727668. If the white 
dwarf actually has a mass as high as is indicated by the 
Doppler boosting technique (i.e., 0.27 M©) then it is dif¬ 
ficult to understand how the white dwarf can still be so 
thermally bloated. If, on the other hand, we take the un¬ 
certainty seriously, and allow for a mass as low as, say, 
0.22 M©, then there is still a discrepancy between the 
inferred cooling age of « 950 Myr, and the apparent lo¬ 
cation of the A star in the HR diagram which is still near 
the ZAMS. 

The growing statistical ensemble of these systems 
should make it worthwhile to (I) make renewed attempts 
to measure the RVs in four of the remaining bloated hot 
white dwarf systems in order to better determine M^d, 
and (2) directly model the binary evolution of such sys¬ 
tems with a code that can evolve two stars simultane¬ 
ously, e.g., MESA (Paxton et al. 2011). 

Note Added In Manuscript: After this work was com¬ 
pleted we became aware of an impressive collection of 
18 similar bloated hot white dwarf systems found in 
the WASP photometric database (Maxted et al. 2013; 
2014; and references therein). Two of these have accu¬ 
rate measurements of the mass, radius, and Teff for the 
bloated hot white dwarf and show multi-periodic pulsa¬ 
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